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Abstract 



We study the unitary time evolution of the entropy of entanglement of a 
one-dimensional system between the degrees of freedom in an interval of length 
t and its complement, starting from a pure state which is not an eigenstate 
of the hamiltonian. We use path integral methods of quantum field theory 
as well as explicit computations for the transverse Ising spin chain. In both 
cases, there is a maximum speed v of propagation of signals. In general the 
entanglement entropy increases linearly with time t up to t = £/2v, after which 
it saturates at a value proportional to £, the coefficient depending on the initial 
state. This behavior may be understood as a consequence of causality. 
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I. INTRODUCTION 



Understanding the degree of entanglement of systems with many degrees of freedom, 
such as quantum spin chains, is currently one of the most challenging problems connecting 
quantum information science and statistical physics. In the last few years several measures of 
entanglement have been proposed [1,2], and have been calculated (analytically in the simplest 
instances, otherwise numerically) in the ground state of several many-body systems. Despite 
a wide investigation, there is still no consensus on the correct measure of entanglement in the 
ground state of a many-body system. In this paper we will consider only the entanglement 
entropy, since it is the measure most readily suited to analytic investigations, especially if 
quantum field theory methods are employed [3]. 

The entanglement entropy is defined as follows. Suppose the whole system is in a pure 
quantum state |\&), with density matrix p = and an observer A measures only a 

subset A of a complete set of commuting observables, while another observer B may measure 
the remainder. A's reduced density matrix is pa = Tr B p. The entanglement entropy is just 
the von Neumann entropy Sa = — Tr^ pa log p\ associated with this reduced density matrix. 
It is easy to check that Sa — Sb- For an unentangled product state, Sa = 0. Conversely, Sa 
should be a maximum for a maximally entangled state. One of the most striking features 
of the entanglement entropy is the universal behavior displayed at and close to a quantum 
phase transition. In fact, close to a quantum critical point, where the correlation length £ is 
much larger than the lattice spacing a, the long-distance behavior of the correlations in the 
ground state of a quantum spin chain are effectively described by a 1+1 dimensional quantum 
field theory. At the critical point, where £ diverges, the field theory, if Lorentz-invariant, is 
also a conformal field theory (CFT). In the latter case, it has been shown that, if A is an 
interval of length £ in an infinite chain, Sa ~ (c/3) \og£ [4,3], where c is the central charge of 
the corresponding CFT. In higher dimensions, the critical entanglement entropy Sa scales 
like the area separating A and B [5,2,3]. These results have been verified by analytic and 
numerical calculations on integrable quantum spin chains and in more complicated systems 
[2,6-10]. When the correlation length £ is large, but finite, it has been shown that for 
1-dimensional systems Sa = -A(c/6)log£ [3], where A is the number of boundary points 
between A and B (the 1-d "area"). This has been confirmed in some integrable models 
with A — 1 [3] and A = 2 [7,8]. Finally the dependence of the entanglement entropy 
on the central charge has been the starting point for understanding the connection between 
irreversibility along renormalization group trajectories (Zamolodchikov's c-theorem [11]) and 
loss of entanglement [2,12]. 

A less investigated situation is how (one measure of) entanglement evolves when the 
system is prepared in a state that is not an eigenstate of the corresponding hamiltonian. 
Only a few works considered some entanglement measures in rather specific models [13,10], 
thus universal features as in the ground state have not yet emerged. In this paper we mainly 
consider a scenario that is simply realized experimentally, for example in ultracold atoms 
[14,15]. Consider a quantum system at T = with an hamiltonian H(X) depending on 
a tunable experimental parameter A (for simplicity one may think of the Ising chain in a 
transverse magnetic field A [16], but the scenario is more general). The system is prepared 
in a pure state \ipo), which corresponds to an eigenstate of H(X ) with Ao ^ A. Then, for 
example, at time t = the parameter is suddenly quenched from A to A. In general we 
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expect that j^o) is not an eigenstate of H(X), and thus the system evolves unitarily according 
to the equations of motion given by H(X), in the absence of any dissipation (as we assume 
throughout the paper). In general, the system does not relax to its ground state, as it is the 
case in the presence of dissipation. Within this scenario we consider the time-evolution of 
entanglement entropy. 

In this paper we concentrate on the case of one space dimension, using two comple- 
mentary methods which have been successfully applied to the time-independent case. The 
first uses the path integral representation and the methods of quantum field theory. It 
is straightforward formally to write a path integral representation for the reduced density 
matrix 

PA (t) = Tr B e- iH ^ )^o\e iHWt (1.1) 

and then, following [3], to write the moments TtaPa^) as path integrals over an n-sheeted 
surface with branch cuts along A. To proceed further however, we assume that H (A) corre- 
sponds to a continuum conformal field theory. As in any continuum field theory, the energy 
fluctuations in the state \ip ) are then divergent. The simplest way of regularizing these 
while keeping the calculation tractable is to modify the state according to the prescription 
IV'o) - > e~ eH ^\ipo), which filters out the high energy components. 

With this modification, we find that when A is an interval of length £, in the limit when 
£/e,t/e ^> 1, the entanglement entropy Sa grows linearly with t, up to a time t ~ £/2 (in 
units where the speed of the elementary excitations is unity). Thereafter it saturates at a 
value Sa ~ nc£/12e. These results, and the crossover time t ~ £/2, can be understood in 
terms of causality applied to entangled left- and right-moving pairs of quasiparticles emitted 
from the initial state. 

As a complement to the general conformal field theory calculations, we calculate the 
time-dependent entanglement entropy Sa for the simplest integrable model: The Ising chain 
in a transverse magnetic field. The reason for this calculation is twofold. On the one hand, it 
is a practical check of the CFT results just summarized and can help understanding problems 
connected with the regularization e. On the other, it allows exact calculations even in the 
non-critical (massive) phase, thus generalizing the CFT results. 

We consider the time evolution from the ground state at a given transverse magnetic 
field h , according to the equation of motion given by H{h) with h ^ h . We obtain the 
remarkable analytic results (equation (3.19)) that for t — > oo, Sa is proportional to £ for 
each pair of h and ho (as long as h ^ ho). This confirms the CFT results and further shows 
that the linear dependence on £ for large times is not a peculiarity of the critical evolution. 
Thus for each given h, there is an arbitrary large entanglement entropy in the asymptotic 
state, contrarily to the ground-state case, where infinite entropy is allowed only at criticality. 
Furthermore we obtain that lim^oo Sa/£ is bounded from above by 2 log 2 — 1 (the value 
found for ho = oo, h = 1). This maximum value of the entanglement entropy is less than 
£log2, a rigorous bound for a generic state of £ spins one-half. An interesting (and not 
understood) feature of the leading term of Sa for t ^> £ — > oo is that it is invariant under 
the exchange h <-> ho- 

For finite times we find a linear increase of Sa with t, for all h ^ ho, not only at criticality. 
The crossover between the linear and non-linear regime happens always at t* — £/2, for large 
£. However for t > £/2, Sa does not immediately saturate to the asymptotic value, but is 
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a slowly increasing function of the time. In Sec. IV, we ascribe this difference between 
lattice and CFT results to the non-linear dispersion relation of quasiparticles emitted from 
the initial state. When both h and ho get closer to the critical field, the time dependence 
becomes sharper and more similar to the CFT results. 

The layout of this paper is as follows. In the next section, we discuss the time dependence 
on the entanglement entropy within CFT, and in Sec. Ill we consider the quantum Ising 
model in a transverse field (some cumbersome details of the calculation are postponed to 
two appendices). In Sec. IV we show that the results obtained may be explained in terms 
of a physical picture invoking causality. Finally in Sec. V we discuss some points that need 
further investigation. 



II. TIME EVOLUTION OF ENTANGLEMENT ENTROPY IN QUANTUM 

FIELD THEORY 

A. Path integral formulation 

Suppose (in the notation of [3]) we prepare the system in a state \tpo(x)) and unitarily 
evolve it with the hamiltonian H. The matrix elements of the density matrix at time t are 



(r(x")\p(t)W(x')) = Z^(r(x")\e- 



-UH—eH 



\^o(x)){i/} (x)\e 



+itH-eH 



W{x>)) . 



(2.1) 



We have included damping factors e~ eH in such a way as to make the path integral absolutely 
convergent. We shall see at the end of the calculation whether it is justified to remove them. 

Each of the factors may be represented by an analytically continued path integral in 
imaginary time: the first one over fields ip(x,r) which take the boundary values ipo(x) on 
t — — 6 — it and if)"(x) on r = 0, and the second with i/j(z,t) taking the values if)'(x) on 
r = and ipo(x) on r = e — it. This is illustrated in Fig. I. 
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FIG. 1. Space-imaginary time regions for the density matrix in (2.1) 
The normalization factor 
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which guarantees that Tr p = 1 is given by the euclidean path integral over imaginary time 
2e, with initial and final conditions both equal to ipo(x). This is the same as sewing together 
the two edges in Fig. 1 along r = 0. 

As in [3], the reduced density matrix is giving by sewing together along r = only 
those parts of the x-axis corresponding to points in B, leaving open slits along A, and Tr p n A 
is given by sewing together n copies of this in a cyclic fashion. Thus 

Tr fi = Z n /Z?, (2.3) 

where Z n is the path integral on an n-sheeted surface, where the edges of each sheet corre- 
spond to imaginary times —T\ and t 2 , analytically continued to T\ — e + it and t 2 = e — it, 
and the branch points lie along r = at the boundaries points between A and B. Finally, 
as derived in [3], the entanglement entropy is given by 

(2.4) 

n=l 



B. Conformal mapping 

Now consider the case when H is critical and the field theory is a CFT. Under the 
renormalization group any translationally invariant boundary condition is supposed to flow 
into a boundary fixed point, satisfying conformal boundary conditions. Thus we may assume 
that the state |^ ) corresponds to such boundary conditions on sufficiently long length scales. 
However this means that the conclusions will be valid only asymptotically. 

For real r the strip geometry described above may be obtained from the upper half- 
plane by the conformal mapping w = (2e/7r) log z, with the images of the slits lying along 
aigz = 7TTi/2e. The result for Z n /Z™ in the upper half z-plane, for an arbitrary number of 
branch points, was given in [3]. It has the form of a correlation function of local operators 
$± n located at the branch points, which transforms simply under conformal mappings. 

First consider the case when A is a slit of length £ (and B is the rest of the real axis.) In 
that case we can take the image of the branch points to be Z\ = p^e 1 ^/ 26 and z 2 = pe tnT1 ^ 2e 
where p = e 7rf//4e . From [3] in the z-plane we have 

/ 1 - 1 1 - 1 \ 2n ^i 

Tr p\ = (KZ-n) - c n ( 1 j! 1 "! 2 r 2 "! 1 , =t ) > ( 2 - 5 ) 



where A n = (c/24)(l — 1/n 2 ) in the notation of [3]. This formula is supposed to be valid in 
the euclidean region, when all the separations are large. Under the conformal mapping, 

(SnMS-nM) = \w'(z 1 )w'(z 2 )\- 2nA -(^ n (z 1 )^ n (z 2 )) . (2.6) 

After some algebra and continuing to t\ = e + it, we find 

Tr^m , MOA^ ( ^ + ^ + 2 COSh(7rt/£) ^ (?7) 
Trp A (t)~cM2e) ^ _ 2 ) • (2.7) 
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In the case where £/e and t/e are large this simplifies to 



c n (vr/2e) 



4nA„ 



-,7r£/2e _|_ g7rt/e 
piz£/2e . pirt/e 



2nA n 



(2.8) 



Differentiating wrt n to get the entropy, 



S A (t) ~ < 



IT27 ( t>£ / 2) - 



(2.9) 



that is increases linearly until it saturates at t — £/2. The sharp cusp in this asymptotic 
result is rounded over a region \t — £/2\ ~ e. 

However we see that e enters in an essential way. The reason is that, in a continuum field 
theory (as compared with a quantum spin model) a state like j^o) has infinitely large mean 
energy (as well as divergent energy fluctuations). In order to make sense of the result it is 
necessary to filter out the high-energy components of the state. Within the path integral 
approach, this is most easily enforced with a cut-off function e~ eE . To compare with results 
from a lattice spin model we should presumably take e to be of the order of the lattice 
spacing. The linear behavior in t for t < £/2, the break at t — £/2, and the saturation at 
a value oc £ all agree with our exact results for the transverse Ising spin chain, in Sec. Ill, 
although there are other differences in detail. 



III. DYNAMICS OF ENTANGLEMENT ENTROPY IN THE QUANTUM ISING 

CHAIN 

As a complement to the general CFT calculation just presented, in this section we de- 
scribe how analogous results can be found in an analytically tractable model. We consider 
the Ising spin chain in a transverse magnetic field, which has a quantum phase transition 
between a ferromagnetic and a (quantum) paramagnetic phase. 

The model is defined by the hamiltonian 

3 

where <rj' 2 are the Pauli matrix acting on the spin at the site j of an infinite chain. The 
quantum critical point is located at h — 1 [16]. We consider the time evolution from an 
initial state l^o) that is an eigenstate of Hj for a field ho ^ h. This experimentally means 
quenching at t — the magnetic field from h to h. We consider only the case h,h Q > 1. 
The generalizations to the case h, h < 1 and to more general spin chains, such as the XY 
model [16] are straightforward and we will not consider them here. 

The determination of the time-dependent state \ip(t)) = e~ iHl ^ t \ipo) (and consequently 
of the entanglement entropy) proceeds with the Jordan- Wigner transformation in terms of 
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Dirac or Majorana fermionic operators. All the details of these computations can be found 
in the Appendix A. 

The final result is that the time-dependent entanglement entropy for I consecutive spins 
in the chain can be obtained (analogously to the ground state case [2]) from the correlation 
matrix of the Majorana operators 



\m<l 



\m<l / 



(3.2) 



We introduce the matrix Tf through the relation (d m a n ) = 5 mn + iTf mn with 1 < m, n < I. 
It has the form of a block Toeplitz matrix 



n n_i ••• ni. 
ni n ; 



n 



-fi 9 1 
-9-i fi 



(3.3) 



with 



where 



9i 



2vr J 



2n 



d(pe tipl e l ° v (cos <E> - % sin <E» cos 2e tp t) , 



2tt 



fl = 7T di P e ^ sin % sin 2e v t ' 
^ Jo 



(3.4) 



\J {h — cos (p) 2 + sin 2 ~ip , 



(ho — cos tp) 2 + sin 2 ip , 
cos (p — h — i sin ip 



sin <E> 



cos $ 



sin <^(/io — h) 



e e° 



1 — cos (f(h + h ) + hh 



(3.5) 



Calling the eigenvalues of Tf as ±iv m , m = !...£, the entanglement entropy is S 
E m =i ^(^m) where if (a;) is 



rr . . l + x. 1 + x 1-x, 1 — x 
H ( x ) = 77- lo S — s 77— lo S -rr - 



(3.6) 



The derivation of all these results is in the Appendix A. 

Note that the first trivial result is that for t — and for h = h , the correlation matrix 
reduces to that of the ground state with field h , as it should. 
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A. Analytical results for t — > oo 



In the limit t — > oo the matrix 1^ simplifies in such a way that analytical calculations 
are feasible for all h, h . In fact, for t — > oo, the time-dependent sine and cosine appearing 
in Eqs. (3.4) average to zero, and Ff has a simpler form with f t — and 



1 /" 27r e - '^ — /i 

9i = — dcpe- tlpl — a-o - ( X ~ cos V?(^ + ho) + hh ) 
^ Jo e <p e p 

Thus it can be written as a block- Toeplitz matrix with 



(3.7) 



(3.8) 



and 



i%) 



o 



/i 1 - cosip(h + ho) + hh 



e l{p — h 1 — cos(p(h + /i ) + ^^o 







(3.9) 



In this case, following the method used for the ground state in Refs. [6,7] it is possible to 
calculate the leading behavior of Se(oo) for large £. Let us introduce [6,7] 



F t = i\I e - Ff , 
L%) = i\h - IL(<p) , 



log 



e(y,x) 



where Ie is the 11 x 2£ identity matrix. 
The determinant of F^ is 



y + x y + x y — x y — x 



log- 



(3.10) 
(3.11) 

(3.12) 



D e (\) = detr e =(-i) e l[(\ 2 -^), 



(3.13) 



so we can use Cauchy theorem to rewrite the entropy as [7] 

S e = lim -L (f d\e(l + e,X)^-lnD e (X), 
e~+o+ Am Jc dX 



(3.14) 



where the contour C depends on the parameter e and includes the interval [—1, 1], so to 
encircle all the zeros of D^{X). This guarantees that the branch points of e(l + e, A) are 
outside C and e(l + e, A) is analytic there. As e goes to zero the contour approaches the 
interval [—1,1]. 

An old short overview about theorems concerning this type of matrices can be found in 
[17]. In particular, it is known that the large t asymptotic behavior of the determinant of 
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2468 10 2468 10 



h h 

FIG. 2. Leading asymptotic entanglement entropy for unit length of the interval, as function 
of h Eq. (3.19). (a) The case ho = oo (full line) is compared with the numerical results (crosses) 
obtained by diagonalizing Tf for t = 200 and using the exact coefficients Eq. (Bl). (b) As in (a), 
but for several finite values of ho. 



a block Toeplitz matrix Ti (under some hypothesis satisfied in our case t = oo) is given by 
the generalization of the Szego lemma [17] 

logdetT* = i— / <fy>log[deta(y>)] + O(logf) , (3.15) 
^ Jo 

where a(ip) is the generating matrix-function of Tt (in our case U(ip)). Thus 

logdetf, = ^-[ ^log[detl%)]+0(log£), (3.16) 
^ Jo 

so (we omit 0(\og£)) 

dlogdetf,(A) _ £ ^ 2A 

aC P. n (l-cosy(fe+fe )+fc/i ( )) 2 \ 6 - 1 ') 

(h 2 +l-2/icos ip)(h%+l-2h cos ip) 



d\ 2ll L \2 (l-cosip(h+h )+hh„) 2 

J u a Tm~i m ._\/t.2 1 1 kz : 



Inserting this expression in Eq. (3.14) we get 

If d 

S e (oo)= lim <p dXe(l + e,X) — IndetlVA) 

2tt 



If £ f Z7T 2X 

lim <b dX e(l + e, A) — / dip r- ,,,,,,,, , 9 . (3.18) 

-^o+ Am In 2ir L \2 (i-cos •p(h+h )+hh ) 2 v ) 



(h 2 +l-2h cos ip) cos <p) 
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Changing the order of the integrals, using the residue theorem, and the property e(y,x) = 
e(y, —x) we get 1 



£ 1 f 2n f 
SAoo) = lim — / dip <b d\ e(l + e, A) 



2A 



^2 (1— cos ifi(h+ho)+hho) 2 

(h 2 +l-2/i cos ip)(h%+l—2ho cos if) 



^dipH - + M + , (319) 

2vr J y ^/(h 2 + l-2/i cos <p) (7ijj + 1 - 2h cos ^) 1 

Let us discuss this result. First, independently of h and /iq, S^(oo)/£ is non- vanishing 
(as long as h ^ ho). Thus the asymptotic entanglement entropy for the Ising chain is always 
linear in £, not only at the critical point as we already knew from CFT. For h — 1, ho — oo 
it reduces to Se(oo)/£ = 2 log 2 — 1, as it can be understood from a very simple argument 
reported in Appendix B. This is the maximum value that Se(oo)/£ can asymptotically 
reach. In Fig. 2 (a) we depict Se(oo)/£ for a quench from h = oo as a function of h: It 
is a monotonous decreasing function of h. In Fig. 2 (b) we depict Se(oo)/£ as function of 
h, for several ho- For h < ho it is a decreasing function, while for h > ho it increases, and 
obviously it is zero at h — h . A further curious (and not understood) feature is that S^(oo) 
is symmetric under the exchange of h and h . 

Note that in Eq. (3.19) corrections of the order 0(\og£) are neglected. These corrections 
are smaller than the leading term, but can be arbitrarily large in the limit £ — > oo. They 
are in principle obtainable using generalizations of the Fisher-Hartwig conjecture [7], and 
are not expected to be symmetric under h <-> h . 



B. Some numerical results 



Now, we present the time-dependence of the entanglement entropy, as obtained by nu- 
merical diagonalization of Tf. We consider several quenches in the next subsections. 



1 . Quench from ho = oo to h = 1 . 

We start our numerical analysis from the case in which the system is quenched from a 
completely uncorrected state (ho = oo) to the critical point (h — 1). The time evolution 
is governed by the critical hamiltonian, so the results we obtain are expected to be well- 
described by CFT. 

The numerical results, as obtained by diagonalizing numerically Tf, are shown in Fig. 
3 for £ = 20,40,60, 100. For all £ considered, Se(t) exhibits just two time-regimes: a linear 



In a similar way one can also calculate the the moments Trp^. In fact, using the results reported 
in the Appendix, it is straightforward to show that logTrp^ is given by equation (3.19) with H 
replaced by H n = log{[(l - x)/2} n + [(1 + x)/2] n }. It follows that Trp^ oc e - fn( - h ' h °^, where f n 
satisfies f n (h,h ) = f n (ho,h). 
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20 40 60 80 100 



t 

FIG. 3. Entanglement entropy for the quench from ho = oo to h = 1, for various I. The dashed 
lines are the leading asymptotic results for large I, cf. Eq. (3.19). The inset shows the derivative 
with respect to the time of Sioo(*)- 

one for t < t* = 1/2 + O(£ ), and a nonlinear one for t > t*. The crossover between the 
two regimes happens at t* ~ £/2, in agreement with the CFT prediction. Note (inset of 
Fig. 3) that the derivative with respect to the time of £ioo(£) is practically constant for 
t < 50, apart from an expected short-time non- universal behavior (that should disappear in 
the limit £ — > oo). For t < t*, all the Sg(t) fall on the same master curve, independently of 
£. However, the entropy does not saturate to the asymptotic value exactly at t* , as in CFT. 
In fact for large t, St(t) is a slowly increasing function, even in the limit of large £. 

The dashed lines in Fig. 3 are the leading asymptotic results for large I as given by 
Eq. (3.19). The actual asymptotic value (obtained by numerical diagonalization) is always 
slightly larger than (3.19), showing that the first correction (most probably of the order 
0(\og£)) is positive. 

2. Quench from ho > 1 to h = 1. 

To understand if the disagreement between CFT and the exact result for t > t* can 
be attributed simply to the particular initial condition (h = oo) or to deeper reasons, we 
consider now the case of a generic initial condition (i.e. finite ho), always evolving according 
to the critical hamiltonian. 

The numerical results for various ho (at fixed I = 60) are displayed in Fig. 4. Even in this 
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t 

FIG. 4. Seo(t) for the quench from ho = oo, 5, 2, 1.5, 1.1 to h = 1. The dashed lines are the 
leading asymptotic results for large i Eq. (3.19). 

case the curves show for t < t* a linear regime, in agreement with CFT. However, decreasing 
ho the time dependence for t > t* becomes sharper, i.e. Se(t) reaches the asymptotic value 
much faster than in the case h = oo. Thus the curves look more similar to the CFT result. 

In Fig. 4 the dashed lines are the leading asymptotic results for large I as given by Eq. 
(3.19). Note that decreasing ho, the positivity of the corrections to Eq. (3.19) is more and 
more apparent. In this case, it can be partially attributed to a non-vanishing value for t = 0, 
that results in a finite shift of the entire curve. 

3. Quench from ho = oo to h > 1 . 

In the lattice calculation we are performing we are not restricted to study the critical 
evolution. Thus, to understand how the entanglement entropy is affected by a non-critical 
evolution, we can simply diagonalize T e A for h ^ 1. As a first case, we consider a quench from 
ho = oo to a non critical value h > 1. The numerical results, for fixed i = 100, are shown 
in Fig. 5. The curves appear very similar to each other, and in particular to the critical 
case. There are two time regimes as in the critical case. The crossover from the linear to 
non-linear regime is always at t* ~ £/2, independently of i and h. A quite unexpected fact 
is that the all the curves fall on the same master curve, when rescaled according to 6*100(00) 
(see the inset in Fig. 5). We have no simple explanation for this. 
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t 

FIG. 5. Sioo(t) for the quench from ho = oo to h = 5, 2, 1.5, 1.01, 1. The dashed lines are the 
leading asymptotic results for large I Eq. (3.19). The inset shows the rescaling of the curves, 
according to the asymptotic value Sioo(oo). 

4- Quench from ho > 1 to h > 1. 

Finally we consider the most general quench into the paramagnetic phase. Some numer- 
ical results for I = 100 are shown in Fig. 6. Due to the double variation of the parameters, 
these figures are not easily readable. The most interesting features are: 

• The crossover from the linear to non-linear regime is always at t* ~ £/2. 

• At fixed ho, when the curves are rescaled according to Se(t = oo), they do not fall on 
the same curve as for h = oo. 

• The dashed lines in Fig. 6 are the leading asymptotic results for large I as given by 
Eq. (3.19). From the figures it is not always clear whether Eq. (3.19) is always smaller 
than the actual value (especially for small h ). We checked in all the cases not enough 
clear graphically (i.e. ho = 1.1 with h = 10, 5, 2 and ho = 1.5 with h = 10, 5) that for 
longer non-depicted times the curves finally cross the value given by Eq. (3.19). This 
shows quite firmly that for all h, ho the leading corrections to Eq. (3.19) are positive. 

• The actual asymptotic values for t — > oo for the cases with h and ho interchanged 
are close, but not exactly equal as in Eq. (3.19). We checked that, for larger I, their 
difference is almost constant and so can be imputed to finite (or 0(\og£)) corrections 
that consequently have not the symmetry under the exchange of h and ho as the leading 
term. 
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FIG. 6. 5*100 (i) f° r a general quench from some generic ho to several h. The dashed lines are 
the leading asymptotic results for large i Eq. (3.19). 



IV. PHYSICAL INTERPRETATION 



The qualitative, and many of the quantitative, features of SU(£) found in the last two 
sections may be understood physically as follows. We emphasize that this is not an ab initio 
calculation but rather a simplified model which allows us to explain these observations. The 
initial state \4>o) has a very high energy relative to the ground state of the hamiltonian H(\) 
which governs the subsequent time evolution, and therefore acts as a source of quasiparticle 
excitations. Particles emitted from different points (further apart than the correlation length 
in the initial state) are incoherent, but pairs of particles moving to the left or right from 
a given point are highly entangled. We suppose that the cross-section for producing such 
a pair of particles of momenta (p',p") is f(p',p"), and that, once they separate, they move 
classically. This will of course depends on H(X) and the state \ipo), and in principle is 
calculable, but we made no strong assumptions on its form. If the quasiparticle dispersion 
relation is E = E(p), the classical velocity is v(p) = dE/dp. We assume that there is 
a maximum allowed speed which is taken to be 1, that is \v(p)\ < 1. A quasiparticle of 
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t 

t 



2t>l 

FIG. 7. Space-time picture illustrating how the entanglement between an interval A and the 
rest of the system, due to oppositely moving coherent quasiparticles, grows linearly and then 
saturates. The case where the particles move only along the light cones is shown here for clarity. 

momentum p produced at x is therefore at x + v(p)t at time t, ignoring scattering effects. 

Now consider these quasiparticles as they reach either A or B at time t. The field at 
some point x' G A will be entangled with that at a point x" G B if a pair of entangled 
particles emitted from a point x arrive simultaneously at x' and x" (see Fig. 7). 

The entanglement entropy between x' and x" is proportional to the length of the interval 
in x for which this can be satisfied. Thus the total entanglement entropy is 



s A (t) 



[ dx' f dx" f dx [ f(p',p")dp'dp"5(x' -x-v(p')t)5(x" -x-v(p")t). 

(4.1) 

Now specialize to the case where A is an interval of length t. The total entanglement 
is twice that between A and the real axis to the right of A, which corresponds to taking 
p' < 0, p" > in the above. The integrations over the coordinates then give max ((v(-p') + 
v(p"))t,t), so that 

/0 poo 
dp' / dp"f(p',p")(v(-p') + v{p")) H(£ - (v(-p') + v{p"))t) + 

dp' / dp"f(p',p") H((v(-p') + v(p"))t - £) , (4.2) 

oo JO 
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where H(x) — 1 if x > and zero otherwise. Now since \v(p)\ < 1, the second term cannot 
contribute if t < t* — £/2, so that S^t) is strictly proportional to t. On the other hand as 
t — > oo, the first term is negligible (this assumes that v(p) does not vanish except at isolated 
points), and Sa is asymptotically proportional to £, as found earlier. 

However, unless \v\ = 1 everywhere (as is the case for the conformal field theory cal- 
culation), Sa is not strictly proportional to £ for t > t*. In fact, it is easy to see that 
the asymptotic limit is always approached from below, as found for the Ising spin chain in 
Sec. III. The rate of approach depends on the behavior of f(p',p") in the regions where 
v(-p') + v{p") — > 0. This generally happens at the zone boundary, and, for a non-critical 
quench, also at p' — p" — 0. If we assume that / is non-zero in those regions, we find a 
correction term ~ —£ 3 /t 2 in the limit where t t*. 

The above result may be generalized to the case when A consists of several disjoint 
intervals, and checked against the corresponding conformal invariance result in the case 
where \v(p)\ = 1. This is discussed in Appendix C. In particular one should note that 
»SU(£) is not always non- decreasing. For example, if A consists of a regular array of intervals, 
each of length £, separated by gaps also of length £, the entanglement entropy oscillates in 
a saw-tooth fashion with a period £. 



We have presented results on the time-dependence of the entanglement entropy Sa, 
starting from a translationally invariant pure state |^o) which is not an eigenstate of the 
time evolution operator, from two complementary perspectives. The first was conformal field 
theory, which applies only in the asymptotic regime, at the critical point, to conformally 
invariant initial states, in theories with a purely relativistic dispersion relation. In order 
to regularize these calculations and maintain tractability we were forced to apply a high 
energy cut-off e~ eE to the state. The second was a solvable lattice model, the Ising model 
in a transverse field, in which we were able to perform calculations starting from a variety 
of initial states, considering both critical and non-critical quenches. In these calculations 
the only dimensional microscopic parameter is the lattice spacing. It should be noted that 
the continuum limit of the critical transverse Ising spin chain is known to correspond to a 
conformal field theory. 

The results of these two sets of calculations have a number of important features in 
common: the entanglement entropy increases linearly with time t (after transients die away 
in the lattice case), up to t* = £/2, in units where the maximum propagation speed of 
excitations is taken to be unity. For t ^> t*, Sa/£ saturates at an asymptotic value. 

However, there are also a number of differences. Partly these may be explained by the 
set-up being different: in the CFT calculation the scale of slope of the linear t-dependence, 
as well as the asymptotic value of Sa/£, is set by e, while for the lattice calculation it is set 
by the lattice spacing as well as the parameter ho specifying the initial state. However, in 
the CFT calculation the discontinuity in the slope is rounded only over a time interval 0(e), 
after which Sa reaches its asymptotic value immediately. This means that the ratio 



V. DISCUSSION 




(5.1) 
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is exactly unity (in units where the maximum speed of propagation is 1.) On the other hand, 
the lattice calculations exhibit a slow increase towards the asymptotic value for t > t*, so 
that R is always less than unity. In Sec. IV however, we argue that this can be understood 
within a simple physical picture and the fact that, on the lattice, there are quasiparticle 
excitations which travel with a group velocity less than the maximum allowed value. This 
also accounts for the existence of a sharp feature in Sa(£) even for non-critical quenches to 
h > 1, which should correspond, in the scaling limit, to a relativistic theory with a mass 
gap. In these cases we find the explicit formula 

R _ r^d P 'j-dff(p',fM-p')+v(p'')] < i 

2jl oo dp'f oo dp"f( P ',p") 

There are however other puzzling features of the Ising results, which may however be 
specific to this model, for example the facts that the limiting behavior lim^oo Sa/(- is a 
symmetric function of h and h, and that the detailed time-dependences of all the quenches 
from ho = oo to any value of h > 1 all appear to collapse all on the same curve when 
suitably rescaled. Clearly these deserve more investigation to see whether they persist in 
other models. 

In this paper, we have restricted attention to quenches from a translationally invariant 
initial state |^o)- It is also interesting to consider other setups, for example when initially 
A and B are completely decoupled and in their respective ground states. This is currently 
under investigation. It should also be noted that the causal picture of Sec. (IV) may be 
simply generalized to higher dimensions, with results which are, however, more dependent 
on the precise geometry. 

Finally we should point out that the increase of the entanglement entropy we have found 
for a single interval A in an infinite system is in no way connected with the second law of 
thermodynamics. As shown in Appendix C, in more complicated situations, it may decrease 
or even oscillate. 
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APPENDIX A: COMPUTATION OF S e FOR THE ISING CHAIN 

The diagonalization of Hj is presented in several standard textbooks (see e.g. [16]). Nev- 
ertheless we briefly present it here to make the notations self-explanatory. The hamiltonian 
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for a chain of length N is 



(JV-l)/2 



(Al) 



J=-(JV-l)/2 



and we consider open boundary conditions. We are interested in the limit iV — > oo. The first 
step in diagonalizing the hamiltonian is to introduce a set of spinless fermion annihilation 
and creation operators through the non-local Jordan- Wigner transformations 



ai 



n 

\m<l 



of - iaf 



af + iaf 



(A2) 



\m<l 



where a] stands for the Hermitian conjugate of di. They satisfy the anticommutation rela- 
tions 

{a}, a m } = Si m , {ai, a m } = 0. (A3) 
Let us now introduce the Fourier transformed fermionic operators 



(N-l)/2 



di 



/ iv, ? 



cue 



(A4) 



l=-(N-l)/2 



with — (AT— 1)/2 < k < (N — l)/2. Since this transformation is unitary, the anticommutation 
relations {o^, d p } = 5k p remain valid. The hamiltonian now takes the form 



(N-l)/2 
fc=-(JV-l)/2 



h — cos • 



Ink 



N 



d\dt — — sin 



% 2nk 
n sin — — 
2 N 



d_kdk + d_ k d k 



(A5) 



A final unitary transformation is needed to cast the hamiltonian into diagonal form. This 
is the so-called Bogoliubov transformation 



where = cos 6^/2, = sin 6^/2 and 

tan 



sin 



2wk 
N 



cos ^A; - h 



(A6) 



(AT) 



Again, due to unitarity of the Bogoliubov transformation, the operators {bj,} follow the 
usual anticommutation relations {b\, b p } = 5k p - Finally, the hamiltonian takes a diagonal 
form 

(JV-l)/2 

Hi= e fc (S+S fe - 1/2), (A8) 

fc=-(JV-l)/2 

where 



ek 



h — cos 



27Tk\' 



+ sim 



2nk 
~N~' 



(A9) 
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The thermodynamic limit is obtained by defining ip = 2irk/N and sending N —> oo, giving 



E i = \ Z ej>% ( A1 °) 



* dip 
with 

e 2 = (h - cosy?) 2 + sin 2 tp , (All) 

where we omit a constant term from the hamiltonian. 

For the time evolution, we closely follow Ref. [15], whose results are mainly based on 
some previous works [18,19]. Since the hamiltonian is translational invariant, only fermionic 
states with opposite momentum k and —k are mixed. We may therefore write the two 
component column vectors 

t>k \ i a / d k 



B k = (^7J , and D k = (J f fc J . (A12) 

The Bogoliubov transformation relating B k and D k is expressed as D k = R x (9 k )B k , where 

Rp(a) = cos^ -Hc^sin^ (A13) 

and cr M are the Pauli matrices with \i = x,y, z. (These are used for conciseness of notation 
and should not be confused with the operators representing the 'spins' of the Ising chain.) 

For t < 0, the system is taken to be in the ground state corresponding to the field ho- It is 
natural to define b' fermions as those which diagonalize the hamiltonian with ho- Similarly, 
9 k and B' k are given by analogy with (A7) and (A6) with h replaced by h . The state |^ ) 
is therefore the vacuum of & fermions. This in particular implies 

W^o) = (J o)=^(^ + l), (AM) 
that in terms of B k fermions reads 

(ij \B k Bl\ij ) = R x (9' k - 6 k ) l -{a z + l)Rl(9' k - 9 k ) = 1(1 + a z cos0 fc + ^sin0 fc ), (A15) 
with (f) k = 9' k -0 k . 

The time evolution of the operators now proceeds according to the hamiltonian (A8), so 
that B k (t) = U k (t)B k (0), where 

t4W=e-^=(^ e l t j=R z (-2e k t). (A16) 

The expectation values at any time can be easily evaluated. Let us introduce the operators 

u t — d\ + cbi , and cu~ = a\ - hi , (A17) 
and the two-component vector operator 

Qk{t) = fe-(o) = ^ (I) ^ fe(t) ' (A18) 
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where we defined oof as the Fourier transform of oof . Using the algebra of 577(2) matrices, 
it is straightforward to find [15] 



1 — sin <p k sin 2e k t —e l6k (cos 4> k + i sin 0^ cos 2efct) 

. e -* e fc ( cos (p k — i sin 0^ cos 2e k t) 1 + sin 0^ sin 2e k t 



(A19) 



Transforming back to real space we have 

< w iX> = a7 E^"^ 1 " ^0 fc sin2e fc t) , 
fc 

K^ 7 ) = ^ Zy* "^- 1 - s in ^sin2e fe t) , 

(wfw/) = E e ifc( '- j) e-^ fc (- cos fc + i sin fc cos 2e k t) . (A20) 
fc 

Taking the thermodynamic limit, the sum over k becomes an integral over ip and 

n n sin 

tanv k — > tantta = , 

cos ip — h 

<f>k^% = & v -e v , (A2i) 

and e k — > e v defined by Equation (All). 

As stated in the text, it is useful to introduce two Majorana operators [2] at each site 
of the spin chain according to equation (3.2). Operators a m are Hermitian and obey anti- 
commutation relations, 

^mi {d m , 0"n\ 25 mn . (A22) 

In terms of the Dirac operators introduced by the Jordan- Wigner transformation they read 

o-21-i =a\ + a h a 2l = -i(a] - a t ) , 

a 2 i-i-ia 2 i „ t a 2 i-i + ia 2 i r\no\ 
oi = g ' a i = 2 ' ^ ' 

and so, in terms of oo ± they are 

0-21-1 = ' = -i^i ■ (A24) 

The two point time-dependent expectation values of a are easily obtained in terms of Eqs. 
(A20) as 

(a 2 id 2rn ) = -(wfw~) , 

(0-2l-l0-2m-l) = (ufUm) > 
{02l02m-l) = -i{ujfUJ+) , 

(a 2 i-id 2m ) = -i(wfu-) . (A25) 
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The entanglement entropy for the time-dependent state is now obtained exactly as in the 
ground-state, following closely Ref. [2]. From equation (A25), the time dependent correlation 
matrix (d m a n ) = 5 mn + iT^ n of the Majorana operators a reads 

■ n n_i ••• rw 
pA _ iii n : 

.n^v.! n 

The coefficients gi and in the limit of an infinite chain are easily derived from equations 
(A20) and (A25), obtaining the results reported in the text equation (3.4). 

As shown in Ref. [2], the entropy of the reduced density matrix p for £ adjacent spins can 
be computed from T A , in spite of the non-local character of transformation (3.2). Indeed, 
the symmetry 




(A27) 



implies that the trace over a given spin at the site j is equivalent to the trace over the two 
Majorana operators at the same site. Thus the density matrix p e can be reconstructed from 
the restricted 2£ x 2£ correlation matrix 

(a m a n ) = 5 mn + i(Tf) mn , m, n = 1, • • • , 21, (A28) 

where Tf is the matrix obtained by T A cutting the last N — £ rows and columns. 
An orthogonal matrix V bring Tf into the block-diagonal form 

(A29) 

This means that the reduced density matrix can be written as a product pi — Qi <S> • • ■ <S> Qe, 
where each qi has eigenvalues (1 ± vi)/2 and entropy S(qi) = H(ui), where H(x) is given by 
equation (3.6). Thus the entropy of pi is the sum of entropies of the £ uncorrelated modes, 

i 

S, = Y,H{vi). (A30) 
i=i 

Summarizing: for arbitrary time, and in the thermodynamic limit (N — > oo), the time 
evolution of the entropy St of the Ising model can be obtained by 

• evaluating Eq. (3.4) numerically for I — 0, • • • , £ — 1, 

• diagonalizing Tf to obtain vi, and 

• evaluating Se using Eq. (A30). 



IT, 



-fi 9i 
-9-i fi 



(A26) 



vrfv T 



7 1 



1 
-1 
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APPENDIX B: SOME ANALYTIC RESULTS FOR THE CORRELATION 

MATRIX 



In some cases it is possible to get analytic expression for the coefficients g t and f h that 
help the numerical computation in the general case. For example, in the limit h — > oo and 
for t — > oo we get 



9i 



1 r 27T 

2Wo 



(h — cosy?) 



2h'+ 2 

J_ 

2h 







I > 1 

I = 

I = -I 

I < -2 



(Bl) 



Using these values in equation (3.4) reduces the difficulty of the numerical calculations for 
generic times. 

For the a critical quench (h = oo and h — 1), the infinite time value of gi reduces to 



f 



2tt 



9i = ^ / Me 



-i<p(l+l) 



(B2) 



In this case, the Toeplitz matrix has only the two diagonals closest to the principal with 
non-zero entries ±1/2. This matrix is nothing but the Laplacian on the lattice (minus 1, 
and without two elements at (2£, I) and (1,2£) that for £ — > oo are irrelevant). Thus the 
eigenvalues for large £ are 

(B3) 



cos— m = !,...£, 



and the entropy is 

hm^M 

«->oo £ 



= - i j2mcos^)^lj\yH(cosy) = 2^og2-^j . (B4) 



APPENDIX C: GENERAL RESULT FOR AN ARBITRARY NUMBER OF 

INTERVALS. 



In this appendix we sketch the CFT calculation of 5a (t) for the case when A consists of 
the union of the N intervals (u2j-i,U2j) where 1 < j < N and Uk < Uk+i- Trp^ is given as 
usual by the ratio Z n /Z'i which has the form of a correlation function 

(Yl $nK-i + in) Yl ®- n (u 2j + in)) , (CI) 

3 3 

in a strip of width 2e. This is related to the upper half z-plane by w = (2e/ir) logz, whereby 
the points + %T\ are mapped to = e 7r ( u fc+ i ' r i)/ 2,E . if we choose X^fc^fe = 0> ^ e jacobian 
factors in the transformation are constant, independent of the u/~. Analytically continuing 
to Ti = e + it sets z k = ■ie 7r(nfc ~*)/ 2e , that is, on the imaginary axis, while the image points go 
over into z k = —ie n ( Uk+t ^ 2e . The correlation function (CI) is a generalization of that in (2.5) 
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and is a product of several factors: the first involves differences (z k — zi) raised to the power 
±2nA n according to whether k — I is even or odd; the second is the same with z k replaced 
by z k ; and the third involves factors like (z k — zj). It can be seen that the t-dependence 
cancels between the first and second type of terms. On the other hand 

_ _ e n(u k -t)/2e _|_ e Tr(ui+t)/2e ^ e 7rmax(u fc -t,u ( +t)/2e (C2) 

where we have assumed that u k ± t ^> e. This expression is to be raised to the power ±2A n 
according to whether k — I is odd or even. Finally differentiating with respect to n and 
setting n = 1 we find 

S A (t) ~ S A (oo) + ^ ^(-^"-'"WK -*,«, + *). (C3) 

If TV is finite (or more generally the are bounded) the second term vanishes for sufficiently 
large t. At shorter times, SU(t) exhibits piecewise linear behavior in t with cusps whenever 
2t = u k — ui, at which the slope changes by ±(nc/6e according to whether k — / is even 
or odd. In the case of an infinite number of regular intervals, with u k = k£, k G Z, SU(£) 
exhibits a sawtooth behavior. 
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